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In this paper, a parametrization describing the kinematical state of the universe via cosmographic 
approach is considered, where the minimum input is the assumption of the cosmological principle, 
i.e. the Friedmann-Robertson- Walker metric. A distinguished feature is that the result does not 
depend on any gravity theory and dark energy models. As a result, a series of cosmographic 
parameters (deceleration parameter go, jerk parameter jo and snap parameter sq) are constrained 
from the cosmic observations which include type la supernovae (SN) Union2, the Baryon Acoustic 
Oscillation (BAO), the observational Hubble data (OHD), the high redshift Gamma ray bursts 
(GRBs). By using Markov Chain Monte Carlo (MCMC) method, we find the best fit values of 
cosmographic parameters in la regions: H = 74.299±4^;, go = -0.386±g;^, jo = — 4.925±?:||f 
and s = -26.404^° 96 4 

which are improved remarkably. The values of go and jo are consistent 
with flat ACDM model in la region. But the value of so of flat ACDM model will go beyond the 
la region. 

PACS numbers: 
Keywords: cosmography 



I. INTRODUCTION 



The kinematical approach to describe the status of universe is interesting for its distinguished feature that it does 
not rely on any dynamical gravity theory and dark energy models. Then it becomes crucial for its potential ability 
to distinguish cosmological models when a flood of dark energy models and modified gravity theories are proposed 
to explain the current accelerated expansion of our universe. This late time accelerated expansion of our universe 
was firstly revealed by two teams' observation of type la supernovae 0,0]. In general, via the Taylor expansion of 
the scale factor a(t) in terms of cosmic time t, the dimensionless coefficients qo, jo and sq named deceleration, jerk 
and snap parameters are defined respectively, for the detailed forms please see Eq. ([8j [TUJ) in the following. For 
convenience, they are dubbed as cosmographic parameters. These cosmographic parameters, which current values can 
be determined by cosmic observations, describe the kinematical status of our universe. For example, the present value 
of Hubble parameter H describes the present expansion rate of our universe, and a negative value of on means that our 
universe is undergoing an accelerated expansion. This kind of approach is also called cosmography cosmokinetics 
[Ell, or Friedmannless cosmology B, Hi- Recently, this approach was considered by using SN in Ref. Q, SN+GRBs in 
Ref. [lj3] and SN+OHD+BAO in [ll|, where the current status of our universe can be read. On the other hand, for a 
concrete dark energy model or gravity theory, when the Friedmann equation is arrived the corresponding cosmographic 
parameters can be derived by simple calculation. As a consequence, the corresponding parameter spaces can be fixed 
from cosmographic parameters space without implementing annoying data fitting procedure. However, the reliability 
of the cosmographic approach depends crucially on how the cosmographic parameter space is shrunk, in other words, 
the improvement of the figure of merit (FoM). That is the main motivation of this paper. In general, when more 
cosmic observational data sets are added to constrain model parameter space, the more degeneracies between model 
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parameters will be broken. Also the FoM will be improved. So, to investigate the current status of our universe and 
to improve the FoM, the cosmographic parameters will be determined by more cosmic observations. When the SN 
and GRBs are used as distance indicators, the Hubble parameter Hq and the absolute magnitudes of SN and GRBs 
are treated as notorious parameters and marginalized. That is to say, SN and GRBs can not fix the current value of 
Hubble parameter Hq. That is what the authors have done in Ref. |9l. llOl] where the cosmographic parameters qo, jo 
and so were investigated. However, the cosmographic parameters permeate in a relative larger space. Of course, to 
describe the kinematical status of our universe well, one has to shrink the parameter space efficiently. Fortunately, 
when the Hubble parameter Ho is fixed as done in Ref. fTlj j . the parameter space is pinned down effectively. When 
the snap parameter so is included, high redshift observations should be added. So, in this paper we are going to use 
SN, BAO, GRBs, OHD to investigate the cosmographic approach. When SN data sets are used, the systematic errors 
are included. The BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, so 
it is helpful to break the degeneracies between parameters. The OHD data sets are used to fix the Hubble parameter 
Ho- Higher redshift data ponits are from GRBs where the correlation parameters are calibrated via cosmographic 
approach synchronously. For the detailed description of these data sets, please see the Appendix lAl 

This paper is structured as follows. In Section llj the definition of cosmographic parameters and basic expansions 
with respect to redshift z are presented, where to consider the convergence issue, the map from z £ (0, oo) to 
y = z/(l + z) G (0,1) is adopted. To the expansion truncation problem, we compare the expansions with ACDM 
model in the range of redshift involved in this paper. The relative departure of Hubble parameter from that of ACDM 
model is up to 20% at the redshift z ~ 1.75. The dfference of distance modulus between the expansion of luminosity 
distance and that of ACDM model is less than 1.6. Section HIT1 are the main results of this paper. To obtain these 
results, the cosmic observational data sets from SN la, BAO, OHD and GRBs and MCMC method are used. The 
detailed descriptions are shown in the Appendix [A] The main points of this paper are listed as follows: 1). BAO and 
OHD are used to shrink the model parameter space 1 . 2). The calibration of GRBs and constraint to cosmographic 
parameters are carried out synchronously. In this way the so-called circular problem is removed. We summarize the 
results in Tab. [I] and Fig- [Hand Fig. [31 Section [IV] is a brief conclusion. 



II. COSMOGRAPHIC PARAMETERS 

The minimum input of the cosmographic approach is the assumption of the cosmological principle, i.e. the 
Friedmann-Robertson- Walker (FRW) metric 



ds 2 



-dt 2 + a 2 (t) 



dr 2 



1 - kr 



r 2 (d0 2 + sin 2 a ^ 2 



(1) 



where the parameter k = 1,0,-1 denotes spatial curvature for closed, flat and open geometries respectively. In this 
paper, we only consider the spatially flat case k — 0. 
The Hubble parameter H{z) can be expanded as 



dH 

H( Z )=Ho + - 
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where the subscript '0' denotes the value at the present epoch and z = l/a(t) — 1. Via the relation 

1 
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(i + z)H(zy 
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1 After our work, the papers used BAO and OHD appeared in arXiv: J. Q. Xia, et. al, arXiv: 1103.0378 and S. Capozziello, et. al, 
larXiv:110X3 096 
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one has 
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where the cosmographic parameters are defined as follows 
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Then the Hubble parameter can be rewritten in terms of the cosmographic parameters as 

H{z) = H {1 + (1 + q )z + (jo - ql)z 2 /2 + [3g 3 + 3g 2 - jo (3 + 4g ) - s ]z 3 /6 



(6) 

(7) 

(8) 

(9) 

(10) 
(11) 

(12) 



For a spatially flat FRW universe, the luminosity distance can also be expanded in terms of redshift z with the 
cosmographic parameters 



d L (z) = cH^ 1 {z + (1 - q a )z 2 /2 - (1 - g - 3g 2 + jo) z 3 /6 

+ [2- 2g - I5q 2 - 15g 3 + 5j + 10g o jo + «o] z A /24: + ...} . 
Via the relation d A (z) — di,(z)/(l + z) 2 , one has the expansion of d A (z) 

d A (z) = cH Q 1 {z ~ (3 + q )z 2 /2 + {11 -j + 7qo + 3q 2 )z 3 /6 

+ (-50 + 13j - 46g + 10j <Zo - 39gg - 15g 3 + s )z 4 /24 + ...} . 



(13) 



(14) 



To avoid problems with the convergence of the series for the highest redshift objects, these relations are recast in 
terms of the new variable y = z/(l + z) [H Q 

H(y) = Ho {1 + (1 + go)y + (1 + g + Jo/2 - q 2 /2)y 2 + (6 + 3j + 6g - 4g j - 3gg + 3g 3 - s ) 2/ 3 /6 

+ (l + go-2.7ogo + 3go72- So /2)y 4 + 0(y 5 )} (15) 

d L (y) = Cj ff - 1 {j/ + (3-go)y 2 /2 + (ll-jo-5g + 3g2) 2/ 3 /6 

+ (50-7jo-26go + 10g jo + 21gg-15g 3 + S o)2/ 4 /24 + 0(y 5 )} (16) 

d A (y) = cHq 1 {y - (1 + q Q )y 2 /2 - (1 + Jo - 9b - 3g V/6 

+ (-2 + j + 2g + lOjogo - 3q 2 - 15g 3 + s )y 4 /24 + 0(y 5 )} (17) 
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With this new variable, z £ (0,oo) is mapped into y £ (0,1). And the right behavior for series convergence at 
any distance can be retrieved in principle fl2l [l3j . When the convergence problem is solved, one has to concern 
the expansion truncation issue. Of course, with higher orders expansion, more accurate approximation would be 
obtained. However, in this way, one has to introduce more model parameters beyond Hq, qo, jo and Sq. How to 
keep the balance between the free model parameters (or expansion truncation) and comic observational data points 
is another complicated problem. That is beyond the scope of this paper. But we'd like to point out that the way 
out may be the so-call Bayesian evidence method. In fact, we can show the deviations of the expansions from ACDM 
model. For illustration, with fixed value of fl m o = 0.27, the relative departure of Hubble parameter from ACDM 
model (the left panel) and differences of distance modulus to ACDM model (the right panel) are shown in Fig. [1] 
Actually, in the redshift range (z £ [0, 1.75], please see Tab. |H]) of the observational Hubble parameters, the relative 
departure of ACDM model is up to ~ 20% which is almost the same of order of error bars of OHD. In the right panel 
of Fig. [TJ the difference of distance modulus between the expansion of luminosity distance and that of ACDM model 
is shown. At high redshift y ~ 1, the departure is larger up to 4. In the redshift range of this paper, z £ [0, 9], the 
difference of distance modulus is less than 1.6. So, up to the fourth oder of y, these expansions are safe. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.2 0.4 0.6 0.8 1.0 

y y 



FIG. 1: The Hubble parameter departure and differences of distance modulus from ACDM model, where fl m o = 0.27 is fixed. 

As the reader has noticed the Taylor expansion is up to snap parameter sq, with these cosmographic parameters 
the Hubble parameter is of the order z 3 . However, <1l{z) and <Ia(z) are of the order z 4 . This is really from the fact 
that the Hubble parameter has contained one order derivative of time t. When it is up to the same order of di,(z) and 
cIa(z), an extra new parameter has to be introduced. So we will classify the data sets on hand into two cases with 
(Case I: SN+BAO+GRBs) or without (Case II: SN+BAO+GRBs+OHD) the observational Hubble data. Another 
reason is that the cosmic observational data sets of SN and GRBs do not have constraint to Hubble parameter Ho. 
That can be seen clearly from the left panel of Fig. [5] in this paper. So, to fix the current value of Hubble parameter, 
the OHD data sets should be added. The reader can also see that the BAO data set is helpful to shrink the parameter 
space. 



III. RESULTS AND DISCUSSION 

In our calculations, we have taken the total likelihood function L oc e~ x ' 'I 2 to be the products of the separate 
likelihoods of SN (with systematic errors), BAO, GRBs and OHD. Then we get \ 2 

X 2 = X 2 SN + X 2 BAO + XGRBs + XOHD . (18) 

where the separate likelihoods of SN, BAO, GRBs, OHD and the current observational data sets used in this paper 
are shown in the Appendix [A"l 

In our analysis, we perform a global fitting to determine the cosmographic parameters using the MCMC method. 
Our code is based on the publicly available CosmoMC package [l4[ . The results are shown in Table U and Figure 
[2J And the evolution curves of the Hubble parameter and distance modulus with respect to redshift z are shown in 
Fig. [3] where the best fitted values of model parameters are adopted from the third row of Tab. Q] 

One can clearly see that when the observational Hubble data are used the lex error parameters space is shrunk 
remarkably. Put in other words, the figure of merit is improved tremendously. It is really from the fact that the 
Hubble parameter H is expressed in terms of z or y with combined cosmographic parameters coefficients. Also, from 
the second row of Table [l] one has noticed that the BAO data set is helpful to break the degeneracy and shrink 
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y.dtu_ 0999 
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TABLE I: The results of xLn, H , go, jo and s in Case I (SN+BAO+GRBs) and Case II (SN+BAO+GRBs+OHD), where 
d.o.f denotes the degree of freedom, a and b are parameters from Amati's correlation of GRBs, for their definition please see 
Eq. IAT31 





FIG. 2: The 1-D marginalized distribution and 2-D contours of model parameter spaces with la, 2a regions. Left Panel: Case 
I: SN+BAO+GRBs. Right Panel: Case II: SN+BAO+GRBs+OHD. 
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FIG. 3: The Hubble parameter and distance modulus with respect to redshift z, where the best fitted model parameter values 
in the third row of Tab. [J are adopted. 

the parameter space. We can test the reliability by comparing the result with spatially flat ACDM model. For the 
spatially fiat ACDM model, we can easily find the corresponding deceleration, jerk and snap parameters respectively 
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(19) 
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When fl m Q varies in the range fi m o G [0,1], qa and sq will be in the ranges qo G [—1,0.5] and sq G [—3.5,1] 
respectively. For comparing the best fit values of cosmographic parameters in Case II with the spatially flat ACDM 
model, where the same data sets combination is used to constrain the flat ACDM model, one finds the corresponding 



result: fL, 



0.270 



+0.0403 
-0.0355' 



-9.398; 



0.0708 j , 
-0.0723 dna 



1.602 
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One can clearly see that for the best fit value 
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of £l m o = 0.270 in flat ACDM model the derived qo = —0.595 and jo = 1 are consistent with the results obtained 
from cosmographic approach in ltr region. However, the value of sg = —0.215 of flat ACDM model is out the range 
of the la region of cosmographic approach. As discussed in Ref. [Tl|, once the parameterized deceleration parameter 
q(z) = qo + qiz/(l + z) [15] is known, one can find the relation qi = —qo — 2q^ — jo- Also one can find other interesting 
relations, for example the relations between the modified gravity theory, DGP brane world model, w = constant, 
CPL parameterized equation of state of dark energy [l3| and cosmographic parameters were investigated in Ref. fl6l | , 
see also in Ref. (Tlj . 



IV. CONCLUSION 



In this paper, the cosmographic approach is reconsidered by using cosmic observational data which include SN 
Union2, BAO, GRBs and OHD via MCMC method. We find the best fit values of cosmographic parameters in 
la ranges: H = 74.299t£!^, q = -0.386±g;^|, jo = — 4.925±|J|f and s = -26.404tj^ 4 which are improved 
remarkably. Comparing with the spatially flat ACDM model, one can find out that the derived values of qo and jo in 
flat ACDM are consistent with the results obtained from cosmographic approach in la region. But the value of Sq of 
flat ACDM model is out of the la region of cosmographic best fit value. As investigated, the BAO data set is helpful 
to shrink the parameter space. When the OHD data sets are added, the parameters space is improved remarkably. 
The reason is from the fact that the Hubble parameter H is expressed in terms of z or y with combined cosmographic 
parameters coefficients. In summary, the main points of this paper are that 1). BAO and OHD are are helpful to 
shrink the parameter space. 2). The calibration of GRBs and constraint to cosmographic parameters are carried out 
synchronously. It is away from the so-called circular problem. 
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Appendix A: Cosmic Observational Data Sets 



1. Type la Supernovae 

Recently, SCP (Supernova Cosmology Project) collaboration released their Union2 dataset which consists of 557 
SN la [Tjj. The distance modulus n(z) is defined as 

Hth(z) = 51og 10 [d L (z)] + (j,q, (Al) 

where dz,(z) is the Hubble-free luminosity distance Ho(Il(z) / c = Ho(1a(z)(1 + z) 2 /c, with Hq the Hubble constant, 
and fj,o = 42.38 — 51og 10 h through the re-normalized quantity h as Hq = lOO/i km s~ Mpc~ x . Where di,(z) is defined 
as 



d L {z) = (1 + z)r(z), r{z) = — — ==sinn 



dz' 



(A2) 



/o E(z')_ 

where E 2 (z) = H 2 (z)/Ho- Additionally, the observed distance moduli ^ bs{ z i) of SN la at z, are 

Uobs{zi) = m obs (zi) - M, (A3) 

where M is their absolute magnitudes. 

For the SN la dataset, the best fit values of the parameters p s can be determined by a likelihood analysis, based 
on the calculation of 

2( 1\T'\ - XT {Vobs(zi) ~ Hth(Ps,Zi)} 2 

X {Ps,M) = 2^ 2 

SN a * 

\ - {51og 10 [J L (p s ,z»)] - m bs(zj) + M'} 2 

SN 1 
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where M' = /io + M is a nuisance parameter which includes the absolute magnitude and the parameter h. The 
nuisance parameter M' can be marginalized over analytically [l|J as 



/+oo 
exp 
-OO 



dM', 



resulting to 



with 



2 , B 2 ( C 



(A5) 



A = {51og 10 [d L (p s ,Zi)] - m o6s (zi)} • Cov^ 1 • {51og 1Q [di(p s ,«,-)] _ ^ofesfe-)} 

SN 

B = Cov y 1 ' { 5 log io (P* ' Z J ' )] ~ m ° bs ( Z J } ' 

i 



(A6) 



where Cov^ is the inverse of covariance matrix with or without systematic errors. One can find the details in Ref. 



|17| and the web site 2 where the covariance matrix with or without systematic errors are included. Relation (|A4| 
has a minimum at the nuisance parameter value M' = B/C, which contains information of the values of h and M 
Therefore, one can extract the values of h and M provided the knowledge of one of them. Finally, the expression 



X 2 SN {Ps,B/C) = A-{B 2 /C), 



(A7) 



which coincides to Eq. (|A5|) up to a constant, is often used in the likelihood analysis [TH, [T^]. Thus in this case the 
results will not be affected by a flat M' distribution. It worths noting that the results will be different with or without 
the systematic errors. In this work, all results are obtained with systematic errors. 



2. BAO 

The BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, and measure 
the distance-redshift relation at z = 0.2. BAO in the clustering of the SDSS luminous red galaxies measure the 
distance-redshift relation at z = 0.35. The observed scale of the BAO calculated from these samples and from the 
combined sample are jointly analyzed using estimates of the correlated errors, to constrain the form of the distance 
measure Dy(z) [20h22| | 



D v (z) = 



{l + zfD\{z) 



H(z) 



1/3 



(A8) 



where Da{z) is the proper (not comoving) angular diameter distance which has the following relation with dh(z 

d L (z) 



Da{z) - (i+zr 

Matching the BAO to have the same measured scale at all redshifts then gives [22 

£V(0.35)/£V(0.2) = 1.736 ± 0.065. 

Then, the X 2 BAO (jPs) is given as 

[LV(0-35)/LV(0-2) - 1.736] 2 



Xbao(Ps 



0.065 2 



(A9) 



(A10) 



(All) 



2 http:/ /supernova. lbl.gov/Union/ 
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3. Gamma Ray Bursts 

Following [24j], we consider the well-known Amati's E Pi i—Ei SO correlation [25U28| in GRBs, where E p ^ — Ep i0 b s (l+z) 
is the cosmological rest-frame spectral peak energy, and E iso is the isotropic energy 

E is0 = 47Td 2 L S bo io/{l + z) (A12) 

in which <1l and Sboio are the luminosity distance and the bolometric fluence of the GRBs respectively. Following 
[H|, we rewrite the Amati's relation as 

log^ = a + Mog-i|l- (A13) 
erg 300keV 

In (2t|, the correlation parameters were calibrated via cosmographic approach. Following this method, we take 
correlation parameters a and b as free parameters when GRBs is used as a cosmic constraint. We fit the Amati's 
relation through the minimization x 2 given by [2~H 

N y-a-bx- 

XGRBsiPs) =J2 „2 l h 2„2 *2 ■ ( A14 ) 



where 



= 1 y.i + b2(J l^ + °%s 



= !«S^7 (A15) 



300keV 
E 

J^iso j 

erg 1 + 2 



Vi = log ^ = log + 2 log 4 (A16) 



where d^ is defined as [3C 



rf L =ff (l + z)r(z)/c, (A17) 
and the errors are calculated by using the error propagation law [3ll ] : 



lnlO^i 
7y ' 4 In 105^0,1 



(A18) 
(A19) 



Here N = 109 GRBs data points are taken from [32|]. The x 2 is large and dominated by the systematic errors, and 
the statistical errors on a and b are small. In general the systematic error a sys can be derived by required \ 2 = v (th e 
degrees of freedom) (24|. Here, we take the value of a 1 sys = 0.324 from Table 1. of the case of Q m o = 0.27 in Ref. [33j j . 
In fact, the concrete value does affect the results concluded in this paper. At last, the total error is af o4 = <J 2 stat + a 2 sya . 
It would be noticed that in our case, the best fit value of a will be less than 2 log(c/i?o) m the definition of luminosity 
distance dt = (1 + z)r(z) [30j . 



4. Observational Hubble Data 

The observational Hubble data are based on differential ages of the galaxies |34|. In (35[ , Jimenez et al. obtained an 
independent estimate for the Hubble parameter using the method developed in [34j . and used it to constrain the EOS 
of dark energy. The Hubble parameter depending on the differential ages as a function of redshift z can be written in 
the form of 

H(z) = --^-^ (A20) 
1 + z dt 

So, once dz/dt is known, H(z) is obtained directly (36|. By using the differential ages of passively-evolving galaxies 
from the Gemini Deep Deep Survey (GDDS) [37j and archival data [38l - l43j . Simon et al. obtained H{z) in the range 
of 0.1 < z < 1.8 [36|]. In [44| . Stern et al. used the new data of the differential ages of passively-evolving galaxies at 
0.35 < z < 1 from Keck observations, SPICES survey and VVDS survey. The twelve observational Hubble data from 
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z 





0.1 


0.17 0.27 


0.4 


0.48 


0.88 


0.9 


1.30 1.43 


1.53 


1.75 


H(z) (km s^MpcT 1 ) 


74.2 


69 


83 77 


95 


97 


90 


117 


168 177 


140 


202 


la uncertainty 


±3.6 


±12 


±8 ±14 


±17 


±60 


±40 


±23 


±17 ±18 


±14 


±40 



TABLE II: The observational H[z) data [Mil]. 



[H, Ell!!] are list in Tabic HI Here, we use the value of Hubble constant H = 74.2 ± 3.6km s ^pc 1 , which is 
obtained by observing 240 long-period Cepheids in [45[ . As pointed out in [45| , the systematic uncertainties have been 
greatly reduced by the unprecedented homogeneity in the periods and metallicity of these Cepheids. For all Cepheids, 
the same instrument and filters are used to reduce the systematic uncertainty related to flux calibration. In addition, 
in [46[ , the authors took the B AO scale as a standard ruler in the radial direction, called " Peak Method" , obtaining 
three more additional data: H(z = 0.24) = 79.69 ± 2.32, H(z = 0.34) = 83.8 ± 2.96, and H(z = 0.43) = 86.45 ± 3.27, 
which are model and scale independent. Here, we just consider the statistical errors. 
The best fit values of the model parameters are determined by minimizing 

= t [g " (Ps; 1rf° bs(Zt)]2 > 

where p s denotes the parameters contained in the model, H t h is the predicted value for the Hubble parameter, H b s 
is the observed value, a(zi) is the standard deviation measurement uncertainty, and the summation is over the 15 
observational Hubble data points at redshifts Z{. The OHD was firstly used to constrain cosmological model in [It}. 
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